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ABSTRACT 


c_ 

The  Ris/a  National  Laboratory ,  Roskilde,  Denmark,  atmos¬ 
pheric  puff  dispersion  aodel  has  been  tested  for  an  atmos¬ 
pheric-marine  environment.  This  three-dimensional  model 
simulates  the  release  of  Guassian  pollutant  puffs  and  pre¬ 
dicts  their  concentration  as  they  are  diffused  and  advected 
downwind  by  a  horizontally  homogeneous,  time-dependent  wind. 
Atmospheric  characteristics  such  as  turbulence  intensity, 
potential  temperature  gradient,  buoyant  heat  flux  and  maxi¬ 
mum  mixing  depth  have  been  considered.  Model  predicted  pol¬ 
lutant  concentrations  have  been  compared  to  airborne  sampled 
observations.  The  effect  of  coastal  turbalence  not  observed 
by  t;ne  single  point  meteorological  measurements  made  onboard 
ship  greatly  affects  the  advection  and  diffusion  of  a  plume 
as  it  moves  onshore.  Additional  measurements/predictions 
particular  to  the  coastal  area  will  have  to  be  incorporated 
into  the  model  for  it  to  accurately  predict  the  onshore 
movement  of  pollutants. 
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I.  INTBODOCTION 


The  downwind  transport  and  distribution  of  atmospheric 
pollutants  froa  an  isolated  source  over  land  or  water  has 
becoae  an  iaportant  environments 1  factor  in  today* s  society. 
The  need  to  understand  the  distribution  of  smoke,  unpleasant 
or  potentially  harmful  foreign  gases  and  perhaps  radioactive 
debris  froa  a  nuclear  powerplant  accident  are  becoming  mr 
and  sore  essential  for  industral  operations  and  construct  n 
planning.  The  dispersion  of  such  atmospheric  pollutantf 
coaaonly  modeled  by  a  standard  Guassian  pluae  model  wh.  i 
coaputes  one-hour  average  plume  characteristics. 

The  Beteorology  Section  of  the  Ris0  National  Laboratory, 
Roskilde,  Denmark,  has  recently  developed  a  puff  model  for 
prediction  and  simulation  of  atmospheric:  pollutant  diffu¬ 
sion.  The  model  considers  individual  puffs  of  pollutants 
with  specific  release  rates  that  are  advected  by  a  horizon¬ 
tally  homogeneous  wind  over  a  grid.  The  wind  input  may  be 
either  the  measured  wind  from  a  single  point,  a  spatial 
average  or  a  wind  simulation.  The  model  simulates  the 
instantaneous  plume  characteristics  by  adding  a  group  of 
puffs,  growing  in  size,  as  they  advect  with  the  wind.  A 
Guassian  pluae  model,  on  the  other  hand,  provides  a  time 
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averaged  concentration  pattern  based  on  a  single  tine  aver¬ 
age  wind  vector.  In  the  puff  model,  the  pluae  advects  with 
a  tiae  series  of  actual  wind  data.  Thus,  the  puff  aodel  is 
able  to  predict  tiae  varying  concentration  distributions  in 
actual  changing  wind  conditions,  waking  it  an  appropriate 
tool  for  dynaaical  coaputations  of  downwind  dispersions  of 
pollutants. 

A  basic  comparison  of  a  puff  aodel  siaulation  and  a  typ¬ 
ical  plume  is  illustrated  in  Pig.  1.  Looking  from  above, 
the  instantaneous  behavior  of  a  plume  being  advected  from  a 
source  by  the  wind  is  shown.  The  outer  cone-shaped  contours 
represent  the  outer  limit  of  the  plume  boundary  and  are 
identical  in  both  Figs.  1  (A}  and  (B)  . 

Fig.  1(A)  shows  an  instantaneous  depiction  of  an  actual 
plume.  The  long-term  average  plume  concentration  is  shown 
on  the  extreme  right  as  a  smooth  curve  with  a  maximum  on  the 
central  axis.  Also  shown  is  the  instantaneous  plume  concen¬ 
tration  considered  realistic  but  is  of  such  a  short  time 
scale  that  it  cannot  be  predicted  or  easily  measured. 

The  puff  model  prediction  is  depicted  in  Fig.  1(B).  The 
circles  show  the  boundaries  of  individual  puffs  of  pollu¬ 


tants  released  from  the  source. 


These  puffs  are  advected 


Figure  1.  Instantaneous  Behavior  of  a  Typical  Plume  and  a 
Series  of  Puffs  from  a  Puff  Model 

and  diffused  downwind  by  a  frequently  updated  wind.  The 
long  term  average  concentration  prediction  of  the  puff  model 
is  expected  to  be  identical  to  the  long  term  concentration 
of  Fig.  1  (A) .  The  short  term  average  pollutant  prediction, 
a  Gaussian  curve  shown  on  the  extreme  right,  is  not  com¬ 
pletely  realistic  but  is  a  reasonable  approximation  to  the 
instantaneous  plume  concentration  profile. 
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The  purpose  of  this  thesis  i:  to  evaluate  and  test  the 
general  characteristics  and  capabilities  of  the  Bis0  puff 
■odel.  The  adjustable  input  data  will  be  varied  and  pre¬ 
dicted  results  for  various  input  cosbinations  cospared.  i 
prelininary  conparison  will  also  be  sade  between  predicted 
results  and  data  collected  fron  a  coastal  region  using 
observed  a eteoro logical  forcing  data  as  sodel  input. 


II.  BIS  Of  POP?  lOfiSI, 

A.  GENERAL  CHAB ACTEBISTICS 

The  Bis#  Puff  Model  is  a  three-dimensional  computer 
■odel  used  for  the  prediction  and/or  simulation  of  the  dif¬ 
fusion  and  advection  of  atmospheric  pollutants.  The  puff 
■odel  technique  is  to  simulate  a  pluae  with  Guassian  shaped 
puffs  with  specified  release  rates  within  a  specified  grid. 
The  initial  size  of  the  puffs  is  noraally  one  aeter  in  diaa- 
eter  although  this  can  be  easily  adjusted.  The  amount  of 
material  in  a  puff  is  the  release  rate  tiaes  the  elapsed 
tiae  between  puffs.  Therefore,  a  long  elapsed  tiae  between 
puff  releases  results  in  a  higher  initial  puff  pollutant 
concentration  than  a  short  tiae  interval.  This  should  not 
noraally  be  of  concern  if  an  adequate  balance  is  maintained 
between  grid  size,  advection  speed  and  puff  release  rate. 

The  location  of  the  puffs  on  the  grid  is  determined  by 
computing  their  aoveaent  for  a  finite  time  step  using  a 
■ensured  wind  field.  The  growth  and  buoyancy  of  the  puffs 
are  computed  froa  siaultaneous  specifications  of  atmospheric 
turbulence  intensity  and  stability  and  from  buoyant  heat 
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flux  at  the  source.  ka  inversion  cap  through  which  pollu¬ 
tants  cannot  pass  and  the  source  height  where  pollutants  are 
released  are  variable  and  can  easily  be  adjusted.  Grid  dis¬ 
tances  within  the  model  say  vary  froa  meters  to  kilometers 
and  tiae  durations  froa  seconds  to  hours  are  possible. 

This  puff  aodel  has  the  capability  of  monitoring  a  aaxi- 
aua  cf  twenty-five  sources  of  puffs  and  its  grid  may  contain 
up  to  100  puffs.  k  puff  source  can  be  located  anywhere  on 
the  grid  and  have  a  unique  release  rate,  start  and  stop  of 
release  tiae,  and  heat  production.  (Then  the  center  of  a 
puff  aoves  outside  the  boundaries  of  the  grid  (either  hori¬ 
zontally  or  vertically)  ,  that  particular  puff  is  dropped 
froa  aeaory.  In  this  way  the  aodel  does  not  store  irrele¬ 
vant  puff  inf or nation,  thus  keeping  coaputer  aeaory  require¬ 
ments  to  a  minimum. 

k  variable  to  control  the  amount  of  reflection/absorp¬ 
tion  of  the  pollutant  by  the  surface  is  easily  adjusted  in 
the  puff  model.  Such  a  parameter  is  of  great  value  both  in 
actual  dispersion  problems  and  also  for  gaining  understand¬ 
ings  of  the  pluae/surface  relationship. 

The  aodel  calculates  the  concentration  at  each  grid 
point  by  suaaing  the  contributions  froa  surrounding  puffs 
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for  each  advection  step.  The  grid  concentrations  can 
be  allowed  to  accuaalate  or  siaply  be  updated  with  the  lat¬ 
est  instantaneous  value.  k  ainiaua  grid  concentration  of 
interest  can  be  set  to  reduce  computer  ran  tiae  by  dropping 
concentrations  too  saall  to  be  of  interest. 

The  output  of  the  aodel  contains  periodic  results  of 
puff  locations  and  concentrations  as  well  as  initial  input 
verification.  The  tiae  interval  for  the  periodic  results  is 
adjusted  by  the  input  data.  This  recurrent  lineprinter  out¬ 
put  contains: 

•  X-Y  plane  plots  showing  the  position 
of  the  sources  and  of  puffs  inside  the 
grid, 

•  X-Z  plane  plots  of  puff  positions 
for  evaluating  pluae  rise  for  each 
vertical  level  of  interest,  and 

•  a  table  listing  of  the  grid  point 
concentrations  for  each  level. 

k  coaputer  drawn  contour  chart  of  the  magnitudes  of  the  pol¬ 
lutant  concentrations  is  also  available. 
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When  considering  distance  between  gridpoints  (delta 
X,Y,Z),  only  spatial  resolution  and  coapnter  resources  need 
be  considered.  Calculated  concentration  accuracy  is  not 
related  to  the  grid-point  separation.  To  insure  that  no 
essential  inforaation  on  individual  puffs  is  "hidden" 
between  grid  points,  the  grid  separation  should  be  adjusted 
dependent  upon  the  size  of  the  puffs  at  the  downwind  dis¬ 
tance  of  interest.  Other  specific  aodel  configuration  con¬ 
siderations  are  described  in  the  following  sections.  They 
are  also  discussed  in  more  detail  in  the  model  behavior 
chapter. 

B.  BIHD  FIELD 

Once  a  puff  is  released,  it  is  advected  based  upon  wind 
data  aeasureaents  at  a  single  point  only,  normally  the 
release  point.  This  limits  the  validity  of  the  aodel  to 
situations  where  the  wind  field  and  turbulence  can  be 
assuaed  to  be  horizonally  homogeneous  throughout  the  grid. 
It  is  therefore  important  to  insure  that  the  data  obtained 
from  such  a  single  point  measurement  is  representative  of 
the  wind  structure  for  the  whole  area  of  interest. 

The  wind  data  are  normally  obtained  in  the  fora  of  a 
horizonal  velocity  time  series.  A  vector  sequence  is  formed 
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by  averaging  over  a  convenient  interval.  These  data  are 
read  into  the  model  after  being  segregated  into  turbulence 
classes  as  discussed  in  the  next  section. 

C.  TUBB ULENCE  INTENSITY  AMD  DIFFUSION 

The  growth/diffusion  of  a  puff  depends  upon  the  turbu¬ 
lence  intensity.  To  account  for  this  growth,  the  puff  eodel 
applies  the  theory  for  relative  diffusion  suggested  by  Saith 
and  Hay  (1961). 

The  turbulence  intensity  is  defined  to  be  the  standard 
deviation  (sigea)  of  the  wind  direction  (in  radians) 
sguared.  The  sigea  values  are  collected  for  the  sane  short 
tiae  periods  as  the  wind  speed  measurenents  used  to  advect 
the  puffs.  Therefore,  the  intensity  of  the  turbulence  which 
governs  the  relative  diffusion  of  the  puffs,  can  be  adjusted 
along  with  the  the  advecting  wind  speed  after  each  time 
step,  if  conditions  warrent. 

A  very  low  value  of  turbulence  intensity  (as  0.0002) 
represents  a  small  standard  deviation  (0.9),  normally  a  sta¬ 
ble  atmosphere  and  a  weak  puff  dispersion/diffusion.  As  the 
atmosphere  becomes  sore  unstable,  the  turbulence  intensity 
increases  along  with  an  increase  of  sigma  values  and  plume 
dispersion/diffusion.  While  these  characteristics  are 
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resentative  of  torbolanca  over  land,  they  can  ba  applied  to 
over  water  cases  in  a  broad  sense. 

D.  FLOSS  BISS 

In  the  vertical  direction,  puff-rise  can  be  accounted 
for  by  Briggs  (1970)  pluee  rise  theory.  In  this  case  buoy¬ 
ancy  is  assused  to  be  conserved  (adiabatic  lotion) ,  and 
pressure  forces,  aolecular  viscosity  and  local  density 
changes  are  considered  ssall  and  are  neglected.  The  rate  at 
which  a  puff  rises  as  it  is  advected  downwind  is  a  function 
of  the  buoyancy  flux,  wind  speed,  puff  distance  traveled  and 
stability  of  the  atsosphere.  Pluae  rise  is  considered  sep¬ 
arately  for  each  individual  puff. 

E.  REFLECTION 

The  interaction  of  the  pollutant  with  the  surface  is 
adjustable  and  can  be  easily  changed  in  the  input  data. 
Total  reflection  or  absorption  or  a  fraction  between  the  two 
can  be  used. 

F.  LIMIT  OF  MIXING  DEPTH 

The  effect  of  an  atiospheric  lid  (inversion)  can  be 
applied  in  the  model  to  limit  the  vertical  movement  of  the 
pollutant.  The  model  does  not  permit  the  pluae  to  rise 
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above  this  cap.  then  a  saxiaas  airing  level  is  in  effect, 
it  acts  to  totally  reflect  the  pollutants  in  the  saae  Banner 
as  total  reflection  at  the  surface.  Phis  airing  cap  would 
act  as  an  inversion  when  the  puff  would  be  erpected  to  grow 
auch  aore  readily  in  the  hori zonal  than  in  the  vertical 
direction. 
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in  intensive  field  tracer  study  vas  performed  during  the 
fall  of  1980  and  winter  of  1981  in  the  Santa  Barbara  Channel 
area  of  the  California  coast.  Fhe  work  was  supported  by  the 
Bureau  of  Land  Banageaent  (BLM)  and  perforaed  by  fhe  Envior- 
aental  Physics  Group  of  the  Naval  Postgraduate  School  (NPS) 
and  ieroVironaent  Inc  (AV) ,  Los  Angeles,  Zk.  This  study  was 
designed  to  help  validate  and/or  aodify  Saussian  dispersion 
aodels  for  coastal  use  and  to  build  a  data  base  for  future 
aodel  development.  Air  pollution  models  in  current  use  have 
not  been  adequately  validated  for  the  over-water  regime. 

In  the  experiment,  a  tracer  gas  (SP6)  was  released  sev¬ 
eral  miles  offshore  from  the  NPS  Research  Vessel  (R/V)  Aca- 
nia.  Ambient  gas  concentrations  as  low  as  10  parts  per 
trillion  (PPT)  were  determined  by  an  array  of  land  based 
sensors,  from  a  snail  boat  and  at  various  levels  by  an  air¬ 
craft  equipped  with  a  continuous  SP6  analyser.  A  chart  of 
the  experimental  area  and  locations  of  the  various  platforms 
is  shown  in  Pig*  2. 
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The  aircraft  flew  through  the  pluae  at  various  eleva¬ 
tions  offshore  and  overland.  The  pluae  transect  tracks  per¬ 
tinent  to  this  study  were  Bade  parallel  to  the  coast  approx- 
iaately  one-half  aile  offshore.  The  airborne  saapling 
consisted  of  instaneous  concentrations  (PPT)  at  selected 
points  at  different  levels  over  a  period  of  six  hours.  The 
observations,  recorded  at  locations  10-70  at  altitudes  of  61 
and  91  a  above  HSL  on  January  29,  1980  are  shown  in  Pig.  3. 
Average  concentrations  over  the  noted  tiae  period  are  shown 
at  the  bottoa  of  each  altitude  block. 

The  following  aarine  aeta orological  paraaeters  were 
aeasured  onboard  the  R/V  Acania  while  anchored  approxiaately 
7.4  ka  offshore: 

•relative  wind  speed  *air  teaperature 

•  wind  speed  fluctuation  *dew  point 

•sea  surface  teaperature  *ship  roll 

•sky  cloud  cover  »ship  location 

•  relative  wind  direction 

•  inversion  height  (acoustic  sounder) 

•vertical  teaperature  and  huaidity  profilies 

(shipboard  radiosonde  launch  every  12  hours). 
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DATE:  9/29/80 
ALTITUDE:  300  FEET  MSL 


TIME:  1700-1800  PDT 
1/2  MILE  OFFSHORE 


The  tracer  gas  was  released  at  a  fixed  rate  through  the 
exhaust  of  one  of  the  ships  motor  generator  sets.  The  gen¬ 
erator  was  run  at  a  constant  speed  resulting  in  a  constant 
stack  temperature  and  flow  rate. 

The  above  collected  Meteorological  inforaation  enabled 
the  data  grid  to  be  established  and  the  atmospheric  wind 
field,  turbulence  intensity,  stability  and  buoyancy  flux  to 
be  derived  for  inputs  into  the  puff  model.  The  puff  model 
prediction  based  on  this  actual  data  formed  the  basis  for 
the  model  performance  evaluations  carried  out  here.  From 
this  basis,  different  input  variables  were  adjusted  to  note 
the  effect  on  the  advected  concentrations — both  in  relation 
to  each  other  and  to  the  airborne  measurements. 
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The  Puff  Model  was  run  on  the  MPS  IBM  S/370  Model  3033 
AP  computer  with  two  goals  in  Bind;  (1)  familiarization  with 
■odel  performance  under  actual  conditions  and  (2)  a  coapari- 
son  of  aodel  predictions  with  observed  data.  The  two  goals 
are  interrelated  in  the  sense  that  atmospheric  data  col¬ 
lected  in  the  aforementioned  tracer  study  were  used  to  fora 
an  initial  prediction  of  the  pi uae  dispersion,  and  varia¬ 
tions  of  that  data  were  used  to  evaluate  the  limits  of  the 
aodel.  The  data  used  as  input  to  the  aodel  represented  the 
marine  atmospheric  conditions  as  determined  from  R/T  Acania 
meteorological  data  at  the  time  of  the  experiment. 

Proper  grid  spacing  was  arrived  at  by  considering  puff 
spread,  mean  wind  direction  and  the  geographical  area  of 
interest.  With  the  initial  prediction  in  hand,  data  input 
variables  of  the  aodel  were  adjusted  and  their  effects 
(changes  in  prediction)  noted.  All  aodel  predictions  were 
compared  with  the  aircraft  observed  data. 

A  7. 4  X  4. 3  1cm  downwind  area  of  interest  was  initially 
gridded  into  an  17  x  10  array.  Distances  between  horizonal 
and  vertical  grid  points  were  approximately  435  a  (Pig.  2). 
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Since  many  of  the  aircraft  observation  tiaes  centered 
around  1730  hours  (all  tiaes  are  Pacific  Daylight  riae) , 
aodel  puff  releases  vere  initiated  at  1630.  The  30-ainute 
wind  speed  averages  obtained  froa  data  taken  onboard  the 
ship  between  1630  and  1730  were  4.7  and  4.8  as-*.  The  first 
advected  puffs  would  be  expected  to  arrive  at  the  back  edge 
of  the  grid  slightly  before  1703,  and  by  1730  a  steady  con¬ 
sistent  pluae  would  be  passing  through  the  area  of  the  air¬ 
craft  track.  (The  aodel  showed,  in  fact,  puffs  leaving  the 
back  edge  of  the  grid  slightly  before  30  ainutes  after  puff 
release)  . 

The  average  wind  direction  was  recorded  on  the  ship 
every  15  seconds.  The  standard  deviations  of  the  wind 
direction  (sigaa)  were  coaputed  as  approxiaate  one  minute 
averages.  These,  in  turn,  were  averaged  over  30  minutes  to 
correspond  with  the  30-minute  wind  speed  averages.  The 
sigaa  values  during  the  time  of  interest  were  1.0  and  0.9 
resulting  in  the  very  small  turbulence  intensity  values  of 
.0003  and  .0002. 

A  delta  Z  value  of  33  a  was  used  to  observe  plume  con¬ 
centrations  at  the  altitudes  of  0,  33,  66  and  99  a  above  the 
surface.  These  levels  were  chosed  for  comparison  with  the 
aircraft  transect  altitudes  of  5  1  and  91  o. 
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Pine  scale  vertical  temperature  and  huaidity  plots  were 
drawn  based  upon  radiosonde  soundings  ta)cen  onboard  the 
ship.  The  sounding  taken  at  1735  POT  (Pigure  4)  shows  a 
shallow  unstable  layer  near  the  surface  topped  by  an  inver¬ 
sion  extending  to  near  400  a.  &  80  a  depth  of  the  sizing 
layer  was  subjectively  established.  The  potential  tempera¬ 
ture  gradient  coaputed  by  the  formula 
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3z 


.0098  z 


(1) 


was  found  to  be  1.0  deg  K/100  s. 

Basic  data  to  determine  source  strength  and  heat  emis¬ 
sion  from  the  ships  stack  were  taken  from  Schacher,  et  al 
(1981).  As  previously  mentioned,  SP6  gas  was  released 
through  the  ships  motor  generator  exhaust  at  a  constant 
rate.  The  stack  temperature  was  250  deg  F,  the  flow  rate 
was  7.13  *  10*  a3s"4.  The  SP8  release  rate  was  47.91  lb 
hr**1.  The  top  of  the  ships  stack--considered  to  be  the 
source  elevation—  was  4  m.  The  source  strength  was  con¬ 
verted  to  6.04  gm  s-1  for  input  into  the  lodel. 

Heat  emission  (H)  in  KW  was  determined  by  the  formula 
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H  3  AT  *  ^  *  Cp  *  Plow 


(2) 


where 

delta  T  =  temperature  (stack  -  air)  dag  K 
B  »  dry  air  gas  constant 
*  2. 87* 10*  erg/g  deg  K 
p  s  io3  P(ab)  *  dyne/cm* 

Cp  »  specific  heat  of  dry  air 
3  . 24  cal/g  deg  K 
Flow  *  16.39  *  7.13  *  103  cm3/s 

The  heat  eaission  was  thus  computed  as  15.07  KW. 

Initially  the  model  grid  was  established  after  noting 
the  area  of  maximum  airborne  sampled  concentrations  (between 
points  24  and  43  of  Fig.  3)  and  the  location  of  the  ship. 
It  became  obvious  during  early  model  runs  that,  with  the 
actual  wind  direction  input,  the  model  predicted  plume  was 
being  advected  south  of  the  grid  towards  point  60  on  the 
aircraft  track.  Obviously,  the  steering  wind,  as  measured 
onboard  ship,  was  not  constant  all  the  way  to  the  shore.  A 
northward  turning  of  the  plume  was  indeed  detected  several 
times  during  the  experiments  by  the  aircraft.  To  compensate 
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for  this  effect,  the  source  of  the  plume  release  was  moved 
in  the  aodel  three  grid  spaces  (  1305  a)  to  the  north  so  that 
the  aaziauB  predicted  pluae  concentrations  would  pass 
through  the  areas  of  the  aaziaua  airborne  aeasured  concen¬ 
trations.  Mo  corrections  were  aade  to  the  aodel  predicted 
pluae  concentrations  because  of  this  adjustment.  However, 
one  could  reason  that  the  predicted  concentration  values 
would  be  higher  in  coaparison  with  aeasured  values  since  the 
coastal  turbulence  and  wind  shift — which  would  tend  to  dif¬ 
fuse  the  plume--were  not  considered. 

Initial  aodel  runs  with  the  saall  turbulence  intensity 
classifications  of  .0003  and  .0002  failed  to  show  pluae  con¬ 
centrations  greater  than  1  *  10-**  in  the  grid  at  any  level 
other  than  at  the  source.  Apparently,  the  grid  spacing  was 
too  large  and  the  narrow  pluae  was  advecting  between  the 
grid  points.  In  an  effort  to  locate  the  pluae,  a  combina¬ 
tion  of  model  runs  were  performed  varying  the  turbulence 
intensity  and  grid  spacing  as  shown  in  Table  I. 

In  this  table,  a  airing  level  cap  of  80  a  was  in  effect 
for  the  model  predictions.  Ho  concentrations  above  that 
level  were  allowed  in  the  computations.  As  previously  men¬ 
tioned,  the  aodel  airing  level  cap  totally  reflects  all 
pollutants  back  downward. 
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A  Comparison  of  Predicted  Concentrations  at  the  Surface,  40  and  80  m 
at  the  East  Edge  of  the  Grid  with  Turbulence  Intensities  between  .01 
and  .05.  Mixing  Level  Limit  is  80  m.  E  =  ( *10) .  Numbers  in 
Parenthesis  are  Grid  Numbers  along  Y  Axis  from  Table  II. 
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With  a  grid  spacing  of  435  a  and  a  turbulence  intensity 
of  .05,  a  plume  concentration  covering  four  grid  spaces  at 
the  east  end  of  the  grid  was  produced.  Predicted  pluee  con¬ 
centrations  slowly  decreased  as  the  turbulence  intensity  was 
reduced  to  .02  (atmospheric  stability  increased).  It  an 
intensity  of  .01,  the  concentration  dropped  by  about  five 
orders  of  magnitude.  Normally,  one  would  expect  increased 
concentrations  with  increased  atmospheric  stability.  Per¬ 
haps  the  result  noted  here  is  due  to  the  plume  shrinking 
away  from  a  grid  point  (and  becoming  more  concentrated 
between  the  recorded  grids)  with  the  increase  in  stability. 

To  increase  the  grid  resolution,  the  grid  spacing  was 
reduced  by  half  to  217.5  a  and  again  to  108.75  m.  With  each 
reduction  the  grid  was  reduced  by  half  in  the  MY"  direction 
and  doubled  in  the  "X"  direction  thus  keeping  distances 
between  grid  spaces  equal  in  all  directions.  This  of  course 
greatly  increases  the  computational  requirements.  If  only 
plume  predictions  along  the  back  edge  are  needed  (as  in 
Table  I)  ,  the  downwind  grid  distance  may  be  held  constant  at 
435  m  while  the  horizonal  crosswind  resolution  is  increased. 
In  this  way  many  unnecessary  computations  are  not  made. 
However,  the  increased  horizonal  crosswind  resolution  is 
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computed  over  the  entire  downwind  grid,  which  in  this  case, 
is  not  necessary.  &  aore  satisfying  solution  to  this  prob¬ 
lem  is  to  install  the  capability  of  using  a  variable  resolu¬ 
tion  grid  with  the  aodel  so  that  downwind  areas  of 
particular  interest  can  be  covered  with  a  dense  grid  while 
ether  areas  of  not  so  such  interest  can  be  sparsely  grided. 

In  order  for  the  advection  of  the  pluae  to  remain  on  the 
array  when  increasing  the  horizonal  crosswind  resolution  and 
decreasing  the  area  exposed  on  the  grid,  the  pluae  source 
was  adjusted  along  the  western  boundary  of  the  grid.  The 
relationship  of  the  vertical  grid  points  to  changes  of  the 
source  location  is  shown  in  Table  II.  The  plume  source  for 
each  grid  resolution  is  noted  with  an  arrowhead.  Grid 
points  that  are  aligned  vertically  in  the  table  have  identi¬ 
cal  locations  and  should  have  the  same  predicted  plume  con¬ 
centrations.  is  mentioned  earlier,  an  increase  of  grid 
resolution  does  not  affect  the  predicted  concentration. 
Notice  that  for  the  same  grid  points  in  Table  I,  the  pre¬ 
dicted  concentrations  with  a  turbulence  intensity  of  .01 
remain  constant  with  changes  in  grid  spacing — only  the  grid 
resolution  was  changed. 
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TABLE  II 


The  Relationship  of  th?  I  Axis  along  the  Western  Grid  Edge 
to  changes  of  Grid  Distance  between  435  and  108.75  a. 

DELTA 

I  (a)  Grid  Points  Along  the  I  Axis 

4  35  4  5  6  7  8  9 

217.5  0123  456789 

108.75  0123  456789 

k 

Proa  Table  I  it  is  obvious  that  the  435  a  grid  spacing  is 
too  large  and  that  the  higher  resolution  does  indeed  "see" 
concentrations  that  would  otherwise  be  missed. 

The  problem  of  a  increasingly  narrow  distance  covered  on 
the  grid  as  resolution  is  increased  can  soaetiaes  be  at 
least  partly  corrected  by  reversing  the  1C  and  T  coordinates 
and  adjusting  or  rotating  the  advecting  wind  direction. 
This  can  easily  be  done  with  the  use  of  the  "TURN"  aodel 
input  paraaeter.  This  procedure  soaetiaes  becoaes  necessary 
since  one  of  the  grid  directions  is  United  by  the  width  of 
the  output  printer  paper  to  less  than  or  equal  10  grid 
units. 

The  aodel  input  variables,  neteorological  and  source 
values  were  adjusted  to  note  their  effect  on  plume  concen¬ 
trations.  A  deeper  understanding  of  how  the  model  works  and 
how  the  ataosphere  affects  dispersion  can  also  be  gained  by 


such  adjustaent s.  &  turbulence  intensity  class  of  .05  was 
used,  except  when  studying  intensity  itself,  because  it  had 
previously  demonstrated  a  good  downwind  grid  coverage  of  the 
pluae. 

In  order  to  note  the  effect  of  the  aaxiaua  mixing  level 
on  the  pluae  concentrations,  several  model  predictions  were 
run,  varying  only  the  height  to  which  the  pluae  was  allowed 
to  rise.  Exact  grid  point  reproductions  were  not  possible 
since  the  model  only  allows  the  height  of  the  mixing  level 
to  be  an  integer  aultiple  of  delta  Z.  The  vertical  grid 
spacing  is  therefore  not  equal.  However,  the  antipicated 
trend  of  increased  concentrations  as  the  mixing  level  is 
lowered  is  evident  from  Table  III. 

The  reflection/absorption  of  the  pluae  at  the  surface  is 
controlled  by  the  model  variable  "REFLEX".  Tests  of  the 
extremes  of  total  absorption  (0.0)  and  total  reflection 
(1.0)  were  performed.  The  results  showed  a  50  percent 
reduction  in  pluae  concentrations  at  the  west  end  of  the 
grid  with  total  absorption  compared  to  total  reflection  in 
otherwise  identical  model  runs. 

The  model  has  a  self-imposed  limitation  of  100  puffs 
from  all  sources  on  the  grid.  The  model  will  terminate  if 
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TABLE  III 


Plume  Concentrations  between  Surface  and  99  a  under 
Different  Maximum  flixing  Levels. 


Max  Mixing 


Level 

(5) 

(6) 

None 

.  44E-8 

.  55E-5 

. 43E-8 

. 54E-5 

. 40E-8 

.  50E-5 

. 36E-8 

.  45E-5 

80  m 

. 57E-8 

.72E-5 

. 59E-8 

.  74E-5 

. 58E-8 

.  72E-5 

0 

0 

30  m 

.  64E-8 

.81E-5 

. 65E-8 

. 8 IE- 5 

0 

0 

(7) 

(8) 

20E-5 

. 22E-9 

SfC 

19E-5 

. 21E-9 

33  m 

18E-5 

.  20E-9 

66  m 

16E-5 

.  18E-9 

99  m 

26E-5 

.28E-9 

Sfc 

27E-5 

.  29E-9 

40  m 

26E-5 

. 29E-9 

80  m 

0 

0 

>  80  m 

29E-5 

.32E-9 

Sfc 

29E-5 

.  32E-9 

30  m 

0 

0 

>30  m 

this  nuaber  is  exceeded.  A  balance  must  be  Bade  between  the 
rate  at  which  puffs  are  released  froa  the  source  (TACJ)  and 
the  time  it  takes  the  puffs  to  be  advected  across  the  grid. 
A  release  rate  of  one  puff  every  40  seconds  was  predomi¬ 
nantly  used  during  this  study. 

The  turbulence  intensity  variable  was  varied  to  include 
conditions  that  are  aore  unstable.  As  atiospheric  instabil¬ 
ity  increases,  the  pluae  would  be  expected  to  expand 
whereas,  with  stable  conditions,  the  pluae  should  reaain 
narrow  and  highly  concentrated. 
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With,  the  use  of  UPS  contouring  routines  and  the 
subroutine  ‘’DRAW*1/  a  visual  comparison  of  the  plume  disposi¬ 
tion  and  concentrations  is  available.  Since  the  computed 
pluae  concentration  varies  over  many  orders  of  magnitude, 
the  concentration  values  were  converted  to  integer  numbers 
by  multiplying  by  1.  *  10  *3  and  then  talcing  the  logarithm. 
These  logarithas  are  then  saoothed.  Thus,  a  contour  plot 
representing  order  of  aagnituda  concentrations  was  be  pro¬ 
duced.  As  with  the  aodel  variables  "MAPTIH"  and  "KPLANS" 
which  control  the  frequency  and  vertical  levels  of  printer 
plots,  "DRAB"  can  be  called  to  contour  concentrations  at  any 
tiae  period  and  for  any  level  required. 

Pluae  concentration  distributions  for  turbulence  inten¬ 
sities  of  .05,  .10  and  .25,  all  other  variables  constant, 
are  shown  in  Pigs.  5,  6  and  7.  As  expected,  the  plume 
becomes  wider  and  the  concentration  decreases  as  the  turbu¬ 
lence  intensity/diffusion  increases  and  the  atmosphere 
becomes  less  stable. 

In  Pigs.  5-7,  the  pluae  source  wc.s  located  at  grid  point 
(0,6).  The  wide  pluae  in  the  lower  part  of  the  plot  is  not 
real  but  is  a  function  of  tha  smoothing  routine  spreading 
the  early  puffs  more  than  would  be  expected.  Since  the 
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Figure  5 


Orders  of  Magnitude  of  Plume  Concentration 
Turbulence  Intensity  equal  to  .05. 


with 
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smoothing  routine  would  tend  to  saooth  strong  concentrations 
near  the  source,  the  smoothing  should  be  eliainated  if  the 
priaary  interest  is  near  the  source.  actually  one  would 
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expect  the  puffs  to  behave  as  in  Pig.  3,  froa  Mikkelsen 
(1979),  showing  the  relationship  between  the  puff  size  and 
concentration,  the  rate  of  puff  release  (TkO)  and  the 
advecting  wind  speed  0. 


^concentration 


Figure  8.  Relationship  between  Puff  Size,  Concentration, 

Puff  Release  Rate  (r AO)  and  fcdvec  ng  Wind  Speed 
0  {Mikkelsen,  1979)  . 
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Puffs  would  have  to  travel  the  distance  Xain  before  they 
expand  to  a  size  where  they  effectively  overlap  and  fora  a 
solid  pluae.  Proa  Pigs.  5-7,  one  can  see  that  pluae  concen¬ 
trations  have  increased  with  distance  and  that  a  Xain  has 
been  reached  in  the  middle  to  upper  part  of  the  grid.  It  is 
at  this  point  that  the  puff  aodal  would  be  expected  to  accu¬ 
rately  predict  pluae  concentrations.  If  the  area  of  inter¬ 
est  is  before  the  present  Xain,  the  ralease  rate  of  puffs 
would  need  to  be  increased  so  that  the  successive  puffs 
would  overlap  sooner.  Also  noted  that  as  the  turbulence 
intensity  increases,  the  area  of  aaxiaua  concentration  of 
the  pluae  expands  while  the  central  concentration  decreases. 
This  agrees  with  conservation  of  mass  theorey. 

To  appreciate  the  relative  importance  of  the  source 
strength  and  buoyant  heat  flux,  these  variables  (discussed 
in  Chapter  IV)  were  doubled  separately  and  together  and  the 
concentrations  compared  to  the  concentrations  from  the 
actual  conditions.  Little  or  no  changes  in  concentration 
were  noted  when  the  buoyant  heat  flux  was  doubled  and  source 
strength  remained  the  same.  However,  when  the  source 
strength  was  doubled  and  heat  flux  held  constant,  the  grid 
concentrations  doubled  as  expected.  Thus,  under  existing 
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conditions,  the  source  strength  was  critical  to  the  pre¬ 
dicted  pluae  concentrations  while  the  buoyant  heat  flux, 
within  the  range  tested,  was  not  relevant.  The  vertical 
printer  plots  did  show  an  initial  puff  rise  soon  after 
release  due  to  the  initial  heat  release  but  as  the  puff  rose 
and  expanded,  it  soon  reached  the  aabient  teaperature  and 
leveled  off.  The  buoyant  heat  flux  would  probably  be  aore 
iaportant  when  dealing  with  a  saaller  scale  grid  or  greater 
heat  release. 


V.  Dm  C0££££I20N 


Ho  atteapt  was  Bade  to  compare  actual  puff  aodel  concen¬ 
tration  predictions  at  exact  grid  points  to  aircraft  obser¬ 
vations  for  the  following  reasons: 

•  The  aircraft  locations  were  approximations — the  exact 
locations  were  not  known.  Large  differences  in  predicted 
concentrations  are  seen  with  small  grid  separations  as 
evidenced  in  Table  II. 

•  As  noted  in  Pig.  2,  the  aircraft  observations  were  taken 
over  a  period  of  time  at  different  levels — while  the  puff 
model  produced  multilevel  instantaneous  predictions. 

•  As  mentioned  earlier,  the  actual  wind  was  not  constant 
between  the  ship  observation  site  and  the  oppsite  side 
of  the  grid  near  shore.  Since  the  model  advects  the 
puffs  based  upon  ship  observed  wind,  the  behavior  of 
actual  plume  would  be  different  from  predicted. 

•  Calabration  procedures  for  the  SP6  continuous  analyzer 
mounted  onboard  the  aircraft  were  not  available  for 
instantaneous  concentrations  greater  than  1010  PPT. 
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Therefore,  a  question  of  actual  levels  of  SF6 
concentration  in  the  higher  ranges  exists. 

The  puff  eodel  predicted  concentrations  are  expressed  in 
g/m3  while  the  aircraft  observations  are  shown  as  the  volume 
of  SFsper  unit  voluee  of  air  in  PPT.  &  conversion  between 
the  predicted  and  observed  concentrations  was  obtained  by 
computing  the  partial  pressure  and  molecular  weight  of  SP6 
at  standard  pressure  and  temperature.  A  conversion  of 
g/m3  a  (.63*10“**)  *  observed  concentration  (PPT) 
was  thus  found. 

Generally,  the  aircraft  sampled  concentrations  (Fig.  2) 
show  values  between  100  and  8000  PPT.  Converting  these 
observed  concentrations  to  predicted  concentration  units 
gives  values  between  .eS^IO-*  and  .50*10~7  g/m*.  These 
observed  concentrations  are  much  smaller  than  the  values 
shown  in  Table  II.  Perhaps  this  difference  could  be 
explained  by  the  fact  that  the  puff  model  advected  the  plume 
toward  the  coast  in  the  same  direction  under  the  same  very 
stable  conditions  as  observed  on  the  ship.  Any  considera¬ 
tion  of  increased  turbulence  and  wind  shifts  near  shore 
would  be  expected  to  reduce  the  actual  plume  concentrations 
toward  the  observed  concentration  levels. 
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Increasing  the  turbulence  intensity  to  0.25  and  keeping 
all  ether  variables  constant,  the  concentration  values  would 
decrease  to  the  order  of  aagnitude  of  10~*--closer  to  the 
observed  concentrations.  (This  would  reguire  the  wind 
direction  standard  deviation  to  increase  froa  1  to  28). 
However  the  increased  instability  would  cause  the  pluae  to 
spread  over  a  such  greater  area  (Pig.  7)  than  observed  by 
the  aircraft. 
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VI.  COy^LUglQ  MS  4fiD  SECQ.1H2HD4 TION 5 

The  puff  modal  has  been  demonstrated  to  be  a  versatile 
working  dispersion  model.  Different  coabinations  of  input 
variables  showed  the  expected  reasonable  results.  Ihe  dif¬ 
ferences  between  aodel  predicted  and  aircraft  observed  pluae 
concentrations  do  not  seen  to  be  the  fault  of  the  aodel  but 
aainly  that  of  the  highly  variable  aeteoro logical  conditions 
found  along  a  coast. 

Probably  the  aost  obvious  conclusion  reached  froa  this 
study  is  that  predicting  the  behavior  of  a  pluae  moving  over 
a  aarine  environment  onto  a  coastal  regia n  has  significant 
problems.  In  all  probability,  atmospheric  boundary  layer 
conditions  7.4  kn  offshore  can  be  very  different  froa  those 
observed  in  the  more  turbulent  coastal  region.  The  single 
point  meteorological  measurement  at  the  source  should  not  be 
expected  to  adeguately  represent  plume  characteristics  as  it 
nears  a  meteorologically  variable  coastline.  Additional 
observations  (primarily  wind  speed  and  direction),  or  other 
aeans  of  predicting  the  coastal  meteorological  conditions, 
would  have  to  be  incorporated  into  the  puff  model  to  ade¬ 
quately  handle  this  problem. 
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The  advantage  of  incorporating  variable  grid  spacing 
within  the  puff  model  and  the  obvious  beaifits  have  already 
been  discussed. 

Presently ,  the  sizing  cap  of  the  puff  sodel  is  required 
to  be  located  at  an  integer  Multiple  of  delta  Z.  More  flex¬ 
ibility  in  this  parameter  to  include  any  level,  regardless 
of  delta  Z,  would  be  benifical. 

Along  with  the  puff  locations  shown  on  the  lineprinter 
output,  a  maximum  concentration  level  of  each  puff  would  be 
helpful. 

In  future  experiments,  several  aircraft  tracks  should  be 
made  further  out  from  the  coast  in  an  attempt  to  avoid  the 
turbulent  coastal  region.  Observations  thus  obtained  in  a 
noncoastal  environment  would  help  to  verify  the  model  pre¬ 
dictions  without  the  coastal  influence. 


49 


APPENDIX  A 


MAJOR  SECTIONS  OP  THE  POPP  MODEL 

The  Ris0  Puff  Model  has  been  described  by  (Mikkelsen, 
1979).  The  code  also  is  well  documented  with  comment  state¬ 
ments.  With  that  information  and  the  outline  to  be  provided 
in  this  and  the  following  appendices,  the  computational  and 
input/oatput  procedures  will  be  obvious. 

The  program  and  input  data  are  stojed  on  cards  for  the 
sake  of  permanency.  For  efficient  operational  execution, 
the  program  and  input  data  cards  are  read  on  a  disk  within 
the  computer.  The  model  can  then  be  run  at  will  without 
reference  to  the  original  data  cards.  Minor  changes  can 
easily  be  made  directly  on  the  disk  both  to  the  model  and/or 
data  before  each  execution. 

The  model  can  be  separated  into  the  following  main  sec¬ 
tions: 

A.  Input  Data 

B.  Initial 

C.  Calculating 

D.  Output 

E.  Error  Diagnostics 

P.  Subroutines 
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These  vill  be  described  separately  in  the  following  sec¬ 
tions. 

A.  IHP0T  DATA  S ECHOS 

The  input  data  includes  the  variables  shown  in  Table  IT. 

TABLE  IT 

Input  Data  Tariables  for  the  Puff  Dispersion  Model. 

Sind  History  Potential  Teaperature  Gradient 

Turbulence  Intensity  Buoyant  Heat  Flux 

Grid  Dimensions  Minimum  Concentration  of  Interest 

Mixing  Depth  Reflection  at  Ground  Level 

Source  Locations,  Start/Stop  Time,  Strength,  Heat  Emission 
Number  of  Seconds  between  Advection  Steps 
Number  of  Seconds  between  Printouts/Plots 
Number  of  Seconds  between  Puff  Releases 


The  wind  field  and  stability  class  for  the  current  time  step 
are  read  at  the  start  of  the  calculation  section. 

The  variables  listed  above  are  printed  as  a  input  data 
check  and  a  permanent  record  to  accompany  the  actual  output. 
In  most  cases  the  print  command  can  be  overridden  by  YES/NO 
options. 

B.  INITIAL  SECTION 

Based  upon  the  input  data  from  section  (A),  the  initial 
section  specifies  and  initializes  parameters  to  be  used  in 
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the  calculating  section  and  is  passed  only  once  during  exe¬ 
cution  of  the  model.  The  grid  and  so»e  counters  are  ini¬ 
tialized.  Constants  relating  to  reflectance,  mixing  depth 
and  stability  as  well  as  those  controlling  the  size  of  some 
of  the  loops  within  the  model  are  established.  Parameters 
such  as  number  of  puff  releases  per  second,  number  of  advec- 
tion  steps  per  second  and  number  of  advection  steps  per  puff 
release  are  determined. 

C.  CALCULATI ON  SECTXOH 

Using  current  wind  and  stability  class  data  read  at  the 
start  of  the  calculation  section,  the  model  advects  the  puff 
centers  and  calculates  the  growth  rate  and  plume  rise  of  the 
puffs.  It  removes  the  puffs  that  have  left  the  grid  (hori- 
zonally  and/or  vertically)  .  The  predicted  concentration  is 
computed  at  the  grid  points  to  include  pollutants  from  all 
nearby  puffs. 

D.  OUTPUT  SECTION 

For  time  intervals  designated  by  the  input  data,  printer 
plots  of  the  X-I  and  Y-Z  grid  are  produced.  A  maximum  mix¬ 
ing  level  is  marked  on  the  Y-Z  grid  if  in  effect. 
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These  plots  include  the  source  location  and  a  trace  of  the 
pluae  from  the  release  tiae  to  the  maptime.  Also  printed  at 
this  interval  is  a  X-T  table  of  grid  concentrations  for  each 
vertical  level  of  interest.  These  concentrations  can  be 
either  accuaulated  or  actual  concentrations  at  the  plot 
tiae. 

Added  to  the  puff  aodel  is  a  versatic  plotter  routine  to 
saooth  and  contour  the  grid  aagnitude  concentrations  of  the 
above  tables. 

E.  ERROR  DIAGNOSTIC  SECTION 

If  the  model  is  directed  by  the  input  data  beyond  the 
limits  of  the  design  of  the  program,  the  program  is  termi¬ 
nated  by  way  of  the  error  diagnostic  section.  It  prints 
comments  relating  to  commonly  made  input  errors  enabling  the 
user  to  isolate  problems. 

F.  SUBROUTINES 

The  subroutine  "Sigris"  calculates  the  puff  size  in  the 
horizonal  and  vertical  directions.  It  also  estimates  plume 
rise  associated  with  pollutant  buoyancy. 

The  subroutines  "Ispace”  and  "R space”  are  used  in  the 
framework  of  the  printer  plots. 
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The  subroutine  "Draw"  converts  the  plume  concentrations 
to  a  logarithmic  values,  saoothes  and  then  contours  thea 
using  UPS  inhouse  contour  subroutines.  The  values  are  con¬ 
verted  to  their  logaritha  values  so  that  the  problem  of  con¬ 
touring  over  many  orders  of  magnitude  is  simplified. 
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APPENDIX  B 


PUFF  MODEL  FLOW  CHART 


INPUT 

SECTION 


READ  -  data  input  variables  (primary,  source, 

intensity,  current  stability) ;  averaging 
_  time  step  (tav)  for  wind  data 


INITIAL  SECTION  [  Initialize  constants  j 


READ  &  STORE  all  wind  speed  and  direction 
for  current  intensity  class 


CALCULATION 


SECTION 


yes 


next  wind  vector 


I  RECALL  f ran  storage  one  wind  vector  for  next  tav  | 


next  advection  step 


Release  puff,  execute  advection  step 
_ calculate  concentration 


3 


Z  'advection 
^  steps 
^.GE.  tav. 


-no 


Z  advection 
steps  '>■ 
.EQ.  MAPTIM 


.yes. 


OUTPUT 

SECTION 

Puirf  positions, 
grid  concentra- | 
tions,  contour 


use  of  data 
|  Termination! 
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APPENDIX  C 


POPP  MODEL  CONTRACTIONS 


CH EMIN — Hiniaua  grid  concentration  of  interest 
DELX ,DELT, DELZ--Dist ance  in  aeters  between  grid  points 
DOSE — Allows  the  concentration  aatrix  to  accnaulate 
DTDZ — Potential  teaperature  gradient  (K/M)  (.GE.  0) 

HEAT — Individual  source  heat  eaission  (KWt 
ICOLS — Nuaber  of  coluans  in  grid  (-LE.  10) 

INST — Instantaneous  concentration  aatrix 

ITIHE — Start  tiae 

JROHS — Nuaber  of  rows  in  grid 

KPLANS — Nuaber  of  vertical  levels  in  grid  (includes  surface) 
HA PTIH— Nuaber  of  seconds  between  printer  plots 
NRELSR — Nuaber  of  seconds  to  stop  of  release 
NRHOlT —  Nuaber  of  sources  (.LE.  25) 

NTADV— -Integer  nuaber  of  seconds  between  advection  steps 
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SB FL SC — Reflection  at  ground  level  (0.  none; 1.0  total) 

SOORRR — Number  to  identify  source 

SODRST — Strength  for  individual  source  (ga/s) 

STOPBL — Individual  source  stop  tiae  (s) 

STRTBL-- Individual  source  start  tine  (s) 

TAD — Integer  number  of  seconds  between  puff  releases 
TORN — Angle  of  rotation  of  wind  direction 
IS ODBCS — X  coordinate  of  source  in  grid  units 
YSOUBCE- -T  coordinate  of  source  in  gria  units 
ZH--Liaited  mixing  depth  (a) 
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